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Abstract. 

Significant advances have been made in recent years in the exploration of clustering 
in light nuclei. This progress has arisen not only from the investigation of new sys- 
tems, but also through the development and application of novel probes. This paper 
will briefly review selected topics concerning halo and molecular states in light nuclei 
through examples provided by the neutron-rich Be isotopes. 



INTRODUCTION 



Clustering within nuclei is a widespread phenomenon! which takes on many guises 
across the nuclear landscape. Until relatively recently cluster studies have been 
confined to on or near the line of beta stability where the role of a-clustering has 
long been established [1] . As clustering is expected to manifest itself most strongly 
near thresholds [2] , exotic structures might be expected to form in very neutron (or 
proton) rich systems. Over the last decade the exploration of clustering in nuclei 
far from stability has become technically feasible as demonstrated most clearly by 
the discovery and subsequent probing of the nuclear halo [3]. Whilst an excess of 
neutrons (or protons) may naively be expected to dilute any underlying a-cluster 
structures, theoretical [4,5] and experimental work [6,7] indicate that molecular 
type structures such as a-chains "bound" by valence neutrons may exist. 

In the present paper a number of selected topics concerning the study of halo 
and molecular states in light, neutron- rich nuclei will be reviewed. As they exhibit 
many of the facets of clustering and structural evolution far from stability the 
neutron-rich Be isotopes have been chosen as examples. In parallel, the techniques 
which have been developed to aid in probing the structure of such nuclei far from 
stability will also be discussed. 
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FIGURE 1. Left: Position spectra for ^°Be ions from the p(^^Be,^"Be)d reaction on a (CH2)„ 
target. Right: Angular distributions for transfer to the ^°Be ground and 2+ states. Comparison 
is made with various core coupUng descriptions for ^^Be [10,11]. 

REACTION SPECTROSCOPY 

Single-Nucleon Transfer Reactions 

One of the issues of prime importance in the study of nuclei far from stabihty is 
the extraction of rehable spectroscopic information. For nuclei on or near stability 
light-ion (p, d, ...) induced single-nucleon transfer reactions have long been one 
of the tools of choice. The application of such reactions to nuclei far from stabil- 
ity presents a number of experimental challenges [8,9]. Most significantly, whilst 
the cross sections are moderately high (~1-10 mb/sr) the short halfiives of the 
nuclei of interest dictate that beams of these nuclei be used in inverse kinematics 
reactions. The consequent constraints on the final excitation energy and angular 
resolution limit the target thickness to the order of 1 mg/cm^ and beams of at least 
some 10^ pps are required in order to obtain angular distributions in a reasonable 
measurement time. 

A recent study [10,11] of the structure of the single-neutron halo nucleus ^-"^Be 
via the p(^^Be,^'^Be)d reaction at 35 MeV/nucleon provides a prototypical exam- 
ple of some of the features inherent in such reaction studies. Experimentally a 
magnetic spectrometer operated in a dispersion matched mode was employed to 
detect the heavy ejectile (^°Be). In order to remove the contamination arising from 
reactions on the C in the (CH2)„ target an array of large area silicon detectors 
was used to detect the coincident deuterons (figure 1). With a beam intensity of 



~3xl0^ pps angular distributions of reasonable quality were acquired in some 72 
hours of running (figure 1). 

In addition to the experimental features illustrated by this experiment it should 
be stressed that the analysis is not model independant. For example, not only did 
the weakly bound nature of ^^Be and the d have to be taken into account but realis- 
tic wavefunctions had to be employed - here coupling to the strongly deformed ^°Be 
core needed to be properly accounted for. Indeed, the use of the standard separa- 
tion energy approach to derive the radial form factors lead to an artifically high core 
excited state admixture. As detailed in refs [10,11], careful analysis indicates that 
the ground state of ^■'^Be is dominated (85%) by the admixture corresponding to the 
valence neutron occupying the 2si/2 orbital together with a modest contribution 
(15%) from the core excited (2+(8)i^ld5/2) configuration. 

Single-Nucleon Removal Reactions 

Measurements of one-nucleon removal (or "knockout" ) reactions on hght targets 
have recently been proposed as a spectroscopic tool for high-energy radioactive 
beams [12,13]. This approach has arisen from the development of reaction calcu- 
lations for halo nuclei in which the strong absorption limit [14] and core excited 
states are accounted for [13]. More specifically, the cross sections for the population 
of a given state of the core fragment (/^) may be related to spectroscopic factors 
(C^S'(/^, n/j)) using an extended version [13,15] of the spectator-core model [16] 
to calculate the cross section {asp) for the removal of the nucleon {nlj). 

a{i:) = j:C'S{i:,nlj)as,{nlj,S:/f) (1) 

nlj 

The corresponding momentum distributions are derived within the same eikonal 
formalism [17] or in a simpler fashion using the opaque limit of the Serber model 
[18,19]. As noted above, the integrated cross sections for the population of the 
core excited states are directly related to the associated spectroscopic factors. In 
analogy with transfer reactions, the shape of the core fragment momentum distri- 
butions plays the role of the angular distributions in specifying the / of the removed 
nucleon. Experimentally the core fragment states are identified by the de-excitation 
7-rays emitted in-fiight, whilst a high acceptance magnetic spectrograph is used to 
determine the momenta. 

One of the principal virtues offered by high energy nucleon removal is the appli- 
cability very far from stability where beam rates are low. This ability to function 
with intensities as low as Ipps is a consequence of the large cross sections (~10- 
100 mb), coupled with the high beam energies (^40 MeV/nucleon) which allow for 
the use of thick targets (~100 mg/cm^). 

A particularly clear example^ of the application of the technique may be found 
in the measurement of single-neutron removal from ^^Be [22], whereby the only 



^) Studies of other near dripline and halo nuclei may be found in refs [12,20,21,23]. 
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FIGURE 2. Core fragment longitudinal momentum distributions for one-neutron removal reac- 
tions on C. The solid lines correspond to Glauber model calculations [17]. 

bound core states are the ground (J'^=l/2'^) and 320 keV (l/2~) levels in ^^Be. 
In simple terms, the population of the ^^Be ground state provides a measure of 
the 2hu! admixture present in ^^Be. The results of the experiment, which indicate 
only a ~30% admixture of the p^-configuration in ^^Be, provide direct confirmation 
[24,25] of the breakdown in the N=8 shell closure. 

The potential power of the technique is further illustrated by a recent systematic 
investigation of a broad range of light, neutron-rich psd-shell nuclei [17]. The inclu- 
sive longitudinal momentum distributions which were obtained in a single rigidity 
setting are displayed in figure 2 whereby a number of features are immediately 
apparent. Most notably, the crossing of the N=8 shell and N=14 sub-shell clo- 
sures are associated with a marked reduction in the widths of the core momentum 
distributions (viz, ^^'^^B, ^^'^^C, ^^O and 24,25p^_ rjj^g former effect arises from the 
large z/2si/2 admixtures expected in the ground states of the Z=4-6, N=9 isotones 
[23,26-28], which also persists for N=10, as suggested by recent studies of ^^Be 
[29,30] (see below). A narrowing of the momentum distributions may also be ex- 
pected for N=15 and 16, as in a simple shell model picture the valence neutrons 
occupy the i^2si/2 orbital. Such results demonstrate that coupled with a high ac- 
ceptance, broad range spectrograph, high energy single-nucleon removal reactions 
offer a powerful means to survey structural evolution over a wide range of isospin 
in a single experiment. 

Finally, in the context of high energy reactions, it should also be noted, as de- 
scribed by Tostevin and Al-Khalili [15,31], that few-body Glauber model analyses 
of total reaction cross sections can also provide important constraints on the ground 
state wavefunctions of halo-like nuclei. As presented in the contribution to these 
proceedings by Suzuki, such analyses are now becoming more widespread in their 
application. 
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FIGURE 3. Spin-energy systematics for states observed in ^°Be (from [34]). The trajectories 
for the postulated positive and negative parity molecular bands are indicated. 

NUCLEAR MOLECULAR CLUSTERS 

As alluded to in the introduction, the a-particle plays an important role in the 
structure of light a-conjugate (A=4n) nuclei [1]. This is a direct consequence 
of the strongly bound character of the ^He nucleus and the weakness of the a-a 
interaction, as evidenced by the unbound nature of ^Be. The persistence of cluster 
structures for systems lying away from the line of beta-stability is well illustrated, 
as will be discussed here, by the beryllium isotopes, for which the a-a system may 
be regarded as the basis. 

From a theoretical point of view, prescriptions such as the Molecular-Orbital 
Model (MO) [4] or the Two-Centre Shell Model (TCSM) [32], in which valence 
nucleons are added to the single-particle orbits arising from the two-centre poten- 
tial, provide a successful means to describe the properties of these nuclei. Moreover 
these orbits may be viewed as the analogues of the a and vr-orbitals associated with 
the covalent binding of atomic molecules. The development of fully fledged Anti- 
symmeterised Molecular Dynamics calculations (AMD) [5] is of particular interest 
as the N-nucleon system is modelled without any a priori imposition of an under- 
lying cluster structure. Recent calculations, in particular, suggest the existence of 
two-centred structures in the Be, B and C isotopic chains with valence neutron 
density distributions exhibiting the features of molecular orbitals [5]. 

From an experimental perspective, von Oertzen [6] has compiled systematic ev- 
idence for the existence of dimers in ^~^^Be and ^^^^B. In the case of ^Be, for 
example, the presence of a valence neutron results in a bound (Borromean) sys- 
tem, the ground and excited states of which may be understood in terms of a 
three-body a:n:a molecular structure. In particular, the rotational bands based 
on the ground and low-lying states exhibit large deformations consistent with the 
associated molecular configurations. 

In the case of ^°Be, the experimental evidence for molecular configurations is 




FIGURE 4. Spin-energy systematics for ^^Be ^He + ^He (from [7]). 



rather less well documented. Beyond the established 0^, 2^ and Ij" - 4||" states, the 
locations of the J=5 [6] and 6 members of the of the negative parity band, as well as 
the J=4 and 6 members of the positive parity band have been postulated following 
recent studies of the a-^He breakup of -"^^Be* [33,34] . As displayed in figure 3 [34] , the 
spin-energy trajectories for the bands based on the O2 and IJ" states at ~6 MeV are 
consistent with large deformations as expected for molecular-like a:2n:a structures. 
Moreover, the location of the bandheads just below the threshold for a + ®He decay 
is in accordance with the considerations of Ikeda describing the formation of clusters 
[2] . Further support for the postulated molecular states may be found in the recent 
AMD calculations of Kanada-En'yo and collaborators [5], whereby well developed 
Q;:2n:Q; configurations are predicted for the 0^ and If bands. 

Given the existence of such molecular-type structures in ^°Be, the question nat- 
urally arises as to the existance of similar structures even further from stability. In 
this context the dripline nucleus ^^Be has been investigated. In a recent study^ em- 
ploying the inelastic excitation of an energetic (35 MeV/nucleon) secondary beam 
of ^^Be, evidence has been found in the breakup into ^He+^He of rotational states 
(J=4, 6, 8) in the excitation energy range 10-20 MeV [7]. As illustrated in figure 
4, the inferred momenta of inertia (^^/2$>=0.15±0.04 MeV) and bandhead energy 
(10.8±1.8 MeV) of the observed states are consistent with the cluster decay of a 
molecular structure which may be associated with a:4n:Q; configurations. 

Further experimental support for the molecular nature of these states would 
be the observation of large spectroscopic factors for the associated clusters. The 
measurement of partial decay widths represents, however, formidable experimental 
challenges, though the presence of relatively few decay channels for the states in 
question may facilitate the measurements. Of additional interest is the search for 
in-band gamma transitions which should also furnish information on the degree of 



^) The existence of such structures was hinted at in an earher study by Korsheninnikov et al. 
[35]. 



clustering [36]. Such measurements may be possible in the near future through a- 
pickup reactions - such as ^^C(^'^He/°'^^Be*)^Be - carried out in conjunction with 
high efficiency charged particle and gamma arrays. 

HALO STATES 

Perhaps the most extreme form of clustering is that exhibited by halo nuclei, 
whereby one or more nucleons reside on average well beyond the core potential [3] . 
In the present discussion the two-neutron halo nucleus ^^Be will be used as an illus- 
trative example. In particular, the spectroscopy of ^^'^^Be, continuum excitations 
of ^^Be and the spatial configuration of the halo neutrons will be addressed. 

The tool chosen to investigate ^"^Be in the work described here was a kinemat- 
ically complete measurement of the fragments (^^Be and two neutrons) from the 
dissociation of a 35 MeV/nucleon beam of ^^Be on C and Pb targets. Such a 
measurement is now a relatively standard technique and allows the two-ncutron 
removal cross sections, neutron angular distributions and invariant mass spectra to 
be extracted, and the neutron-neutron correlations to be explored. 

The details of the experiment will not be repeated here as they have already been 
described in refs [29,37,38]. It should be stressed, however, that one of the principle 
problems confronting such measurements is the detection with the highest possible 
efficiency of two beam velocity neutrons at very forward angles with small relative 
momenta and minimal cross-talk. As described in refs [37,38], the use of a highly 
granular array arranged in a staggered configuration coupled with off-line cross-talk 
rejection algorithms (based on kinematic conditions) permit such measurements to 
be undertaken. 

Structure and Continuum Excitations 

The results obtained for the two- neutron removal cross sections, cr_2n, the 
single-neutron angular distributions^, r„, and the associated angle integrated 
cross sections, Un, are displayed in Table 1. The average neutron multiplicities 
{Win = <^n/'^-2n): which havc also been deduced, are instructive in terms of the 
mechanisms leading to dissociation [39] . For a light target the reaction is expected 
to proceed via single-neutron removal (absorption or diffraction) followed by the 
in-flight decay of ^^Be. As approximately equal contributions are expected for 
absorption and diffraction [39] the average neutron multiplicity should be 1.5, in 
accordance with that measured. This scenario is also supported by the single- 
neutron angular distribution for the C target which is well reproduced assuming 
passage via a low-lying resonance in ^^Be [29,40]. In the case of a heavy target, 
nuclear and Coulomb dissociation occur. Given that Coulomb dissociation results 



^) The angular distributions were well characterised by a Lorentzian lineshape. 



2000 



1500 - i 




FIGURE 5. i^Bc + neutron relative energy spectrum from the dissociation of ^*Be on C at 
35 MeV/nucleon [42]. The results of a simulation including states in -'^'^Be at 0.5 (triangles) and 
2.0 MeV (crosses), with r=0.5 MeV are also displayed. 
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r„ [MeV/c] 


G 


0.46±0.04 


0.75±0.10 


1.6±0.3 


75±3 


Pb 


2.3±0.4 


4.0±0.3 


1.7±0.2 


77±4 


Pb(EMD) 


1.45±0.40 


2.7±0.4 


1.9±0.6 


87±6 



TABLE 1. Measured cross sections, average neutron multi- 
plicities and neutron distribution momentum widths for the 
dissociation of ^"'Be at 35 MeV/nucleon [29]. 



in a multiplicity of 2, the average multiplicity for dissociation on Pb should be 
between 1.5 and 2, as observed. 

The possible existence of a low-lying state in ^^Be is of considerable importance 
in modelling the structure of ^^Be [41]; in particular in fixing the ^^Be-n interaction. 
Beyond the single-neutron angular distributions, the relative energy for ^^Be + n 
events may be reconstructed from the measured momenta. The preliminary results 
of such an analysis [42] are displayed in figure 5 together with the results of Monte 
Carlo simulations (see below) in which levels at 0.5 and 2.0 MeV above threshold 
are assumed with widths of r=0.5 MeV. The level at 2.0 MeV is that observed in 
multi-nucleon transfer reaction studies [43,44], whilst the former is consistent with 
a preliminary result reported by the MSU group [45]. Given, as discussed above, 
the appearance of ground states dominated by a valence z/2si/2 configuration for 
the neighbouring N=9 isotones ^^B and ^^C [17,21,23], together with the predicted 
continuation of such an inversion for ^^Be [26], it appears that the ground state is 
most probably J'^=l/2"'". Analysis of the ^^Be - n angular correlations should shed 
further light on this conjecture. 

The enhanced cross section for dissociation on the Pb target (Table 1) is indica- 
tive of a large EMD contribution. Assuming that the nuclear-Coulomb interference 
is small, the C target data may be scaled to estimate the nuclear contribution to 



breakup on Pb [29]. Given a root-mean-square radius of 3.2 fm for ^^Be [30], 
a^l^[Ph) = 0.85±0.07 b and, consequently, a^^^{Pb) = 1.45±0.40 b. 

The invariant mass spectra, reconstructed from the measured momenta of the 
beam and fragments (^^Be and two neutrons) from breakup, are displayed in figure 
6a and b for the C and Pb targets. The EMD spectrum (figure 6c) has been 
deduced following subtraction of the estimated nuclear contribution to reactions on 
Pb and exhibits enhanced strength around 1.5 MeV decay energy (Edecay)- Given 
the complex nature of the response function of the present setup, a detailed Monte 
Garlo simulation, including the influence of all nonactive materials, was developed 
based on the GEANT package [29]. The results shown in figure 6 were obtained 
following the descriptions for dissociation on C and Pb outlined above and after 
filtering through the simulation. In the case of the nuclear induced reactions only 
the lowest low-lying state in ^'^^Be {Eq = 0.5 MeV, Pq = 0.5 MeV) was assumed 
to be populated following the diffraction of one of the halo neutrons. The EMD 
was simulated under the assumption that the energy sharing between the ^^Be and 
the two neutrons was governed by 3-body phase space. As shown in figure 6c, 
the observed EMD decay energy spectrum could be reproduced assuming a Breit- 
Wigner lineshape with Eq = 1.8 ±0.1 MeV and Pq = 0.8 ± 0.4 MeV. Furthermore, 
the corresponding simulations of the single-neutron angular distributions were in 
good agreement with those observed [29] . 

Thompson and Zhukov have examined ^^Be within the framework of a 3-body 
model in which the ^^Be core was treated as inert^ [41] and a number of trial 
wavefunctions developed. Based on the binding energy and matter radius of ^''Be, 
together with the known d-wave resonance at 2.0 MeV in ^^Be, two ^^Be wave- 
functions were favoured (both of which required an s-wave state near threshold in 
-^^Be as suggested above): the so-called D4 wavefunction - 86% z/(2si/2)^ and 10% 
1/(14/2)2; and C7 - 29% 1/(251/2)^ and 67% 1/(14/2)^. The EMD decay energy 
spectra calculated from the corresponding El strength functions [41] are compared 
in figure 6 with that of the empirical Breit-Wigner lineshape deduced from the mea- 
surements. The corresponding integrated two-neutron removal cross sections are 
1.05 b (D4) and 0.395 b (C7) [41], compared to the measured value of 1.45±0.40 b. 
Although the strength is predicted to be concentrated at a somewhat lower energy 
than that observed, a large 1/(251/2)^ admixture to the valence neutrons wavefunc- 
tion is favoured. This conclusion is also supported by the total reaction cross section 
measurement of Suzuki et al. [30] reported in these proceedings. 

A microscopic cluster model has also been used to explore ^^'^"^Be [46]. In the 
case of ^^Be an s-wave state is predicted very close to threshold, whilst the energy 
of the d-wave resonance is well reproduced. Significantly, a strong El transition 
[B(El)^1.2e2fm2] centred at E^ecay = 1-5 MeV is predicted in ^''Be, very close to 
the structure observed experimentally. Analysis of the corresponding energy surface 
suggests, however, that this transition is not associated with a true resonance [46]. 
Gonsequently, the enhanced strength observed near threshold in the EMD invariant 



'^^ An inert core precludes, ab initio, the existence of any simple negative parity resonances. 




FIGURE 6. Left: Reconstructed decay energy spectra for the dissociation of ^"'Be on (a) C, (b) 
Pb and (c) deduced for EMD on Pb. Right: Comparison of the EMD decay energy spectrum 
deduced from the measurement (shaded region) with 3-body calculations (see text). 

mass spectrum may be qualified as a nonresonant "soft" dipole excitation. 

In liglit of tlie breakdown in tlie neutron p-sliell closure in ^^Be discussed previ- 
ously, the importance of p-sd cross shell excitations in the ^^Be ground state should 
be explored. In a first instance, the search for core excited states (^^Be*) in the 
dissociation of ^^Be may provide indications for the presence of such configurations, 
whilst theoretically the inclusion of a deformed core [25], in three-body calculations 
should be undertaken. 

Neutron— Neutron Correlations 

One of the most intriguing questions regarding the description of two-neutron 
halo systems is the degree of correlation between the neutrons. Considering only 
the intrinsic momentum distributions, a spatially compact dineutron will be char- 
acterised by Pcore = '2pn, whcrcas for a system with no correlations Pcore = V^Pn- 
Unfortunately no unperturbed experimental measure of the momentum distribu- 
tions is accessible. 

Perhaps the method the most well adapted to probing spatial correlations is in- 
tensity interferometry [47]. In this approach the relative motion of the two outgoing 
neutrons is governed by the neutron-neutron final-state interaction and quantum 
statistics (Fermi-Dirac), both of which are related to the spatial characteristics of 
the source. Providing that any other effects on the neutron momenta may be ne- 



glected or eliminated in the construction of the correlation function, interferometry 
provides a means to probe the spatial correlations of the halo neutrons [37]. 

Importantly, owing to the low momentum content of the halo neutrons, the 
standard approach to constructing correlation functions (applied, for example, in 
the measurement of leki et al. [48]) is no longer valid [38]. In particular, the narrow 
momentum distributions result in strong residual correlations and consequently a 
significant overestimate of the source size. As a result a new iterative technique has 
recently been developed and applied to the dissociation on a Pb target of beams of 
^He, ^^Li and ^"^Be [37]. Assuming that the timescale between the emission of the 
two neutrons was short (r„„ <100 fm/c), as is most probably the case for EMD, 
neutron-neutron RMS separations of 5.9±1.2 (^He), 6.6±1.5 (^^Li), and 5.4±1.0 fm 
(^^Be) have been deduced. These results appear to exclude the presence of any 
significant dineutron configurations. By way of comparison it is interesting to note 
that the RMS proton-neutron distance in the deuteron is some 3.8 fm. Future high 
statistics experiments emplying a well modelled system such as ^He should allow 
the emission timescale and source sizes to be extracted simultaneously from the 
longitudinal and transverse neutron-neutron relative momenta as well as coherent 
analyses of the neutron-neutron and core-neutron correlations. 

CONCLUSIONS 

In this paper some of the facets of nuclear halo and molecular states in light 
nuclei have been reviewed. In particular, various illustrative examples provided by 
the neutron-rich beryllium isotopes have been discussed. In addition, a number of 
related experimental probes have also been presented. Clearly no single technique 
can furnish a complete description of these nuclei and it will require the appli- 
cation of a broad range of theoretical and experimental tools to obtain a better 
understanding of the phenomena outlined here. 

It is a pleasure to thank my colleagues in the Groupe Exotiques at LPC and in 
the DEMON, CHARISSA and E264/281/295 collaborations who have contributed 
to the work described here. 
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